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Abstract. A class of generalized Randall-Sundrum type II (RS) brane-world models with Weyl fluid are confronted with the 
Gold supernovae data set and BBN constraints. We consider three models with different evolutionary history of the Weyl 
fluid, characterized by the parameter a. For a = Q the Weyl curvature of the bulk appears as dark radiation on the brane, 
while for a = 2 and 3 the brane radiates, leaving a Weyl fluid on the brane with energy density decreasing slower than that 
of (dark) matter. In each case the contribution fl^ of the Weyl fluid represents but a few percent of the energy content of the 
Universe. All models fit reasonably well the Gold2006 data. The best fit model for a = is for fl^ = 0.04. In order to obey 
BBN constraints in this model however, the brane had to radiate at earlier times. 
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INTRODUCTION 

The ACDM model according to which our Universe is a Friedmann-Lemaitre-Robertson- Walker (FLRW) space-time 
with flat spatial sections containing approximately 3% baryonic matter, 24% cold dark matter, the rest being given by 
the contribution of a cosmological constant A seems to be in excellent agreement with current observational data. As 
the dark sector (dark matter and dark energy - a generalization of the vacuum expectation energy represented by the 
cosmological constant) remains unknown, alternate gravitational theories have been advanced. 

The string-theory motivated brane-world models contain our observable Universe as a time-evolving 3-dimensional 
brane embedded in a 5-dimensional bulk. Standard model fields act on the brane, but gravity is allowed to leak into 
the fifth dimension [1], where non-standard model fields could also exist (for a review see [2]). The projection of 
the 5-dimensional Einstein equation onto the brane generates an effective Einstein equation with new source terms as 
compared to general relativity [3], [4]. 

Among them, the energy-momentum squared source term modifies early cosmology [5] and becomes important 
during the final stages of gravitational collapse [6]. It behaves as the dominant source term before the Bing Bang 
Nucleosynthesis (BBN). This quadratic source term is proportional to 1//L. The value of the brane tension X is 
constrained by the deviation from the gravitational Newton-law still compatible with nowadays rigurous experiments 
[7]-[8]. The emerging high value of A and the fast decrease of the square of the energy density of matter implies that 
in a cosmological context this source term can be safely ignored at present time. 

The Weyl curvature of the bulk gives rise to a non-local bulk effect on the brane, appearing as fluid on the brane 
(the Weyl fluid). In the simplest case, when the bulk is static, the Weyl fluid is a radiation field (dark radiation). This 
situation represents an equilibrium configuration, without any energy exchange between the brane and the bulk. BBN 
constraints the amounts of the energy density of the dark radiation as —1.02 x 10^'* < iij < 2.62 x 10^^ [9]. (Here 
£ld is the dimensionless dark radiation energy density parameter.) More generic Weyl fluids were also considered [8], 
[10]-[15]. Depending on how this Weyl fluid evolves, its present day amount can be either negligible or not. It is this 
aspect we wish to consider here, based on our previous analysis [16]-[17]. 

Various brane-world models were confronted with supernova data [18]-[21], however in all these models the 
contribution of Weyl fluid was dropped, assuming it was pure dark radiation during all stages all the cosmological 
evolution. In our analysis we keep the contribution of the (non-radiation like) Weyl fluid. We already gave the analytical 
expression in terms of elliptical integrals for the luminosity distance when the Weyl contribution is small [16] and can 
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be considered a perturbation. ' Then we have tested the models with Weyl fluid characterizing a bulk-brane energy 
exchange, by comparing their predictions with the best available supernova data [17]. Here we present additional 
analysis and strenghten our conclusions. 

Our model consists of a spatially flat FLRW brane embedded symmetrically into a 5-dimensional Vaidya-anti de 
Sitter bulk. The latter has a cosmological constant A, black holes with masses m on either sides of the brane and 
radiation. If the radiation is swithced off, m is constant, and the bulk becomes Schwarzschild-anti de Sitter In this 
configuration the Weyl fluid appears as dark radiation on the brane (with energy density m /a^). Any radiation escaping 
from the brane causes m to vary. The ansatz m = niQa"-, with a = 2, 3, comparable with structure formation has been 
recently advanced [15]. 

The brane coupling constant is related to the coupling constant in the bulk and brane tension A as 6k'^ = k'^X. 
The relation between the brane tension and the bulk and the brane cosmological constant A is 2A = k^X 
We introduce the following dimensionless quantities: 
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Q.ij. Here H is the Huble parameter, p is the matter energy density on the brane and 
a is the scale factor. The subscript denotes the present value of the respective quantities. The Friedmann equation 



where Q.tot = + + 



written in these parameters becomes H /Hq = H^ + Hpaya 
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The luminosity distance for the spatially flat FLRW brane becomes 
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HxaQ/a . At present time this gives 
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The above complicated integral has no analytical form in the majority of cases. However for small H^/ and (as noted 
earlier)) vanishing Hx this integral could be given analytically as d^^^ = -|- il^/Z^, where the coefficient 1^ is 

an analytic expression of elementary functions and elliptic integrals of the first and second kind [16], having different 
forms depending on the actual value of the parameter a. 

We have then compared the predictions of the models characterized by various values of a with the Gold 2006 
supernovae data set [23], for the range —0.1 < < 0.1 up to z = 3. 
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FIGURE 1. For a = the global minimum is at Ha = 0.735, £lp = 0.225 implying = 0.04. The ACDM model is at 
0.\ = 0.725 and fi,„ = 0.275. The white area represents a forbidden parameter range. 



' We note that for a wide class of phantom Friedmann cosmologies similar analytical results in terms of elementary and Weierstrass elliptic functions 
for the luminosity distance are available [22]. 
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The data selects among the brane-world models with a = a global minimum at £2a — 0.735, i2p = 0.225 and 
Q-d = 0.04. We note that the value of Hp is in perfect agreement with the WMAP 3-year data [24]. The l-(7 and 2-(7 
confidence levels are shown of Figl in the Ha — Hp plane. The ACDM model is contained in the 2-(J cofidence level 
at i^A — 0.725, i2p = 0.275. The forbidden (white) region appears because the Friedmann equation combined with 
^tot = 1 gives constraints on the allowed range of Q.^ — Q.p [17]. 

The cases a = 2 and a ~ 3 are represented on Fig2 and Fig3, respectively. In these cases the determination of 
the global minimum becomes unreliable, as the l-a contours become much elongated. Instead of a peak we have 
a "trough", which lies aslope in the Q.j^-Q.p plane. Together with the increase in a the l-a contour becomes more 
elongated and turning counter-clock- wise. 

This feature of the 1-cr contours indicates that with increasing a the models become more likely to be compatible 
with the available supernova data, irrelevant of the exact (but small) value of £2^/. This practically means that any small 
amount of Weyl fluid with a =2, 3 is perfectly compatible with supernova data. 

In the investigated models the Weyl fluid is either dark radiation (a = 0) or describes a situation when the brane 
radiates into the bulk, feeding the bulk black holes (a = 2, 3). In the latter cases the energy density of the Weyl fluid 
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decreases slower, than that of matter, therefore a sizeable Weyl fluid contribution nowasays is perfectly compatible 
with BBN constraints. 

For a = the Weyl fluid evolves as radiation, therefore the BBN constraints can be satisfied only with an 
infinitesimal amount of dark radiation nowadays. The preferred value of £1^ = 0.04 would give too much dark radiation 
in the past. However if we assume that a = is only a recent characteristic of the brane, this equilibrium situation being 
preceded by an epoch in which the brane is allowed to radiate, the BBN constraints can be obeyed with a small, but 
non-negligible amount of dark radiation nowadays. The BBN contraint [9] can be satisfied [17], if the brane radiates 
in the interval (z\, z*)- Assuming zi = 3 the constraint for z* with different value of a gives z* > 61 14.20 for a = 1, 
z* > 155.40 for a = 2, z* > 45.08 for a = 3, z* > 24.01 for a =4. Thus the known history of the Universe can be 
explained if the brane radiates during at least some period of the cosmological evolution. At early times the radiation 
of the brane leads to a black hole, which can further grow during structure formation. 

None of the investigated models, compatible with the available supernova data and structure formation, can be 
excluded by present observational accuracy. The differences among the predictions of the models are however 
increasing with redshift, therefore future measurements of very far supemovae will be able to either support or falsify 
these cosmological models. 
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